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diffusion  coefficients  observed  for  ceramics  subjected  to  pressure¬ 
sintering  or  press-forging  operations  with  interdiffusion  coefficients 
observed  in  diffusion  experiments.  A  magnesium  aluminate  spinel  was 
chosen  as  model  material.  Background,  scope  and  progress  made  to  date  is 
reported  in  Appendix  I,  which  represents  the  Ph.D.  thesis  proposal  of 
Mr.  R.  H.  Smoak.  This  study  is  being  continued  with  funding  from  other 
sources. 

In  addition  to  the  studies  on  basic  mechanisms  in  the  densification 
kinetics  of  ceramics,  considerable  attention  was  also  devoted  to  studies 
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kinetics  of  ceramics,  considerable  attention  was  also  devoted  to  studies 
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of  the  mechanical  behavior  of  ceramic  materials .  The  static  and  cyclic 
fatigue  behavior  of  a  polycrystalline  alumina  was  studied,  which  showed 
the  presence  of  a  strong  cyclic  fatigue  mechanism  in  addition  to  the 
well-known  static  fatigue  mechanism.  Also,  the  fracture  of  glass  was 
analyzed  which  showed  a  direct  relation  between  flaw  size  and  mirror 
radius  of  the  fracture  surface.  In  addition,  the  thermal  stress  fracture 
characteristics  of  ceramics  strengthened  by  surface  compression  was 
investigated  experimentally  as  v;ell  as  theoretically.  The  study  clearly 
shewed  the  advantage  of  surface  compression  strengthening  to  circumvent 
the  severe  problem  of  the  generally  low  thermal  stress  resistance  of 
ceramics.  These  three  studies  on  mechanical  behavior  all  resulted  in  a 
publication  in  the  Journal  of  the  American  Ceramic  Society  of  which  25 
reprints  were  forwarded  to  ARO-D  previously. 

Additional  studies  of  the  mechanical  behavior  of  ceramic  materials 
involved  an  analysis  of  fracture  energies.  In  a  paper  submitted  to  the 
Philosophical  Magazine  (reproduced  in  Appendix  II)  it  was  suggested 
that  measured  values  of  fracture  energy  in  part  consist  of  unrecoverable 
clastic  energy,  such  that  disci'epancies  bet^-;een  surface  free  energy  and 
fracture  energy  do  not  have  to  be  attributed  to  energy  dissipative 
processes  such  as  plastic  floiv.* 

Finally  the  fracture  behavior  of  glass  v;as  analyzed  in  terms  of 
dislocation  mechanisms  v;hich  resulted  in  a  paper  submitted  for  publication 
to  the  American  Ceramic  Society  (reproduced  in  Appendix  III).* 


♦(Note:  Since  Appendices  II  and  III  have  also  been  sxibmitted  in 
manuscript  form  co  AROD,  They  (along  with  App.I)  are 
reproduced  in  their  entirety  in  only  25  copies  of  this 
final  report.) 


Appendix  I 


PRESSIRE  SINTERING  KINETICS  AND  DIFFUSION  IN 
POLYCRYSTALLINE  MAGNESIUM  ALUMINATE  SPINEL 

The  fabrication  of  ceramic  bodies  by  pressure  sintering  makes  it 
possible  to  obtain  high  density,  fine  grained  materials  with  marked  improve^ 
ments  in  mechanical  properties  at  a  temperature  significantly  lower  than  the 
temperature  required  for  a  conventional  sintering  process.  Since  density 
and  microstructure  have  a  large  influence  upon  the  mechanical  behavior  of 
ceramic  bodies  a  knowledge  of  means  to  control  these  factors  becomes  extremely 
important.  In  order  to  improve  the  pressure  sintering  process  it  is  neces¬ 
sary  that  a  better  understanding  of  the  basic  mechanisms  of  densification 
be  developed. 

Densification  Models 

The  model  on  which  most  analyses  of  diffusion-ccxitrolled  densifica¬ 
tion  phenomena  are  based  is  that  of  Nabarro  (1)  and  Herring  (2).  A  diffu- 
sional  mechanism  is  assumed  with  vacancy  migration  occuring  from  regions 
under  tensile  stresses  (vacancy  sources)  to  those  subject  to  compressive 
stresses  (vacancy  sinks).  The  pressure-temperature  regime  where  this  model 
is  valid  will  yield  a  linear  relationship  between  log  e  and  log  with 

a  slope  of  unity.*  This  model,  originally  derived  to  explain  steady-state 
creep  behavior,  has  been  modified  to  account  for  the  rate  of  densification 
during  the  final  stages  of  pressure  sintering. 

*  €  is  the  strain  rate  and  a  is  the  effective  stress  on  the  specimen. 
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Direct  evidence  that  the  Nabarro- Herring  mechanism  of  stress- 
directed  vacancy  flow  is  operative  during  creep  was  obtained  by  Tien  and 
Gamble  (3)  in  their  studies  on  the  influence  of  applied  stress  and  stress 
sense  on  grain  boundary  precipitate  morphology  during  creep.  Before  creep 
studies  were  carried  out  the  micros  true  turn  of  their  tost  specimens*  was 
examined  and  found  to  consist  of  a  regular  array  of  grains  on  the  order  of 
20  ••  50  p,  m  in  size  and  containing  a  uniform  distribution  of  coherent  pre¬ 
cipitates.  After  both  tensile  and  compressive  creep  to  a  point  of  5% 
strain  it  was  found  that  the  boundaries  in  tension  were  denuded  of  preci¬ 
pitates  while  the  volume  fraction  of  precipitates  at  boundaries  in  compres¬ 
sion  was  increased.  By  use  of  the  electron  mlcroprcbe  to  analyze  the  chemistry 
of  the  grain  boundary  areas  it  was  found  that  the  denuded  boundaries  were 
rich  in  chromium  while  the  chromium  content  of  the  boundaries  containing 
high  percentages  of  precipitates  was  depleted.  Since  chromium  is  known  to 
stabilize  the  matrix  at  the  expense  of  the  precipitates  it  was  concluded 
(3)  that  the  Cr  diffused  by  way  of  stress-directed  diffusional  flow  during 
creep  in  the  manner  described  by  Nabarro  (1). 

The  final  stage  pressure  sintering  models  as  derived  from  Nabarro- 
Herring  creep  theory  contain  a  density  function,  f(D)  ,  to  account  for 
the  difference  between  the  effective  stress  and  the  applied  stress  on  the 
specimen.  The  generalized  equation  for  the  final  stage  of  pressure  sin¬ 
tering  is  €  =  K  .  (Here  o^,  the  effective  stress  on  the  specimen, 

*  A  nickel  based  superalloy  of  composition  (wt.  percent)  Ni  -  l6Cr  -  5A1  - 
4Ta. 


is  equal  to  a  f(D)  where  a  is  the  applied  stress;  K  is  a  constant.) 

The  strain  rate  during  pressure  sintering  can  be  equated  to  the  densifica- 

•  -1  dt  1  '’P 

tion  rate  through  the  relation  e  =  "T"  T“  =  ~  where  .*•  and  ■f- 

p.s.  i  dt  P  ■’t 

are  the  instantaneous  relative  density  of  the  compact  and  the  compact 
length,  respectively  f-*). 

In  order  to  determine  a  value  for  the  constant  K  (which  includes 
the  diffusion  coefficient)  from  experimental  data,  it  is  necessary  to  have 
an  expression  for  f  (D).  Various  investigators  (4-7)  have  proposed  such 
expressions  and  demonstrated  their  applicability  with  original  data.  How¬ 
ever  these  expressions  (given  in  Table  I)  have  been  analyzed  by  Coble  (8) 
and  found  to  yield  the  same  result  within  a  factor  of  four.  The  values  of 
the  functions  converge  to  one  as  the  relative  density  approaches  unity  so 
that  at  high  density  the  choice  of  a  density  function  should  not  materially 
affect  the  value  of  the  effective  stress. 

Coble  (9)  has  also  derived  a  pressure  sintering  model  for  the  case 
where  diffusion  via  grain  boundaries  plays  the  major  role  in  the  densifica- 
tion  process.  This  model,  along  with  the  model  based  on  lattice  diffusion, 
is  given  in  Table  II. 

Some  studies  of  creep  of  ceramic  materials  have  been  success¬ 

fully  modeled  by  the  Nabarro-Herr ing  theory.  Coble  (9)  has  noted  that  work 
on  polycrystalline  AI2O2  by  Folv;eiler  (10);  by  Warshaw  (11)  and  by  Beauchamp, 
et  al.  (12)  indicates  that  it  exhibits  typical  Nabarro-Herring  creep  behavior. 
This  was  also  true  for  BeO  as  studied  by  Chang  (13)  and  by  Chandler,  et  al. 
(14).  The  diffusion  coefficients  Calculated  from  the  Nabarro-Herring  model 


HTTSO-o  n  ■<  o  a  C)  o  «• 
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Pressure 


p.s. 

L 


a 


b 


TABLE  II 


Sintering  Models  Based  on  the  Nabarro-Herrlng  Creep  Mechanism 


For  a  lattice  diffusion  controlled  process  (8)  ; 


p.s, 


40D  O 
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For  a  grain  boundary  diffusion  controlled  process  (8) ; 


e 

P 


.s« 


47.5D 

o 

d^T 
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=  steady-state  strain  rate  during  pressure  sintering, 
=  lattice  diffusion  coefficient, 

*  vacancy  volume, 

»  grain  size  (diameter), 

=  applied  stress, 

=  surface  energy  of  the  pore, 

=  pore  radius, 

=  relative  density  of  the  compact, 

=  grain  boundary  diffusion  coefficient, 

=  effective  grain  boundary  width, 

»  Boltzman's  constant,  and 
=  temperature  in  degrees  Kelvin. 
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were  the  same  as  the  values  obtained  from  diffusion  measurements  within  the 
estimated  error  limits  for  their  experiments. 

However,  there  is  a  great  body  of  data  in  the  literature  in  which 
the  diffusion  coefficient  calculated  fay  the  Nabarrc- Herring  model  is  several 
orders  of  magnitude  larger  than  the  intrinsic  diffusion  coefficients.  Sc«i.e 
of  this  information  is  shown  in  Figures  1  and  2  which  were  taken  from  an 
article  by  Vasilos  and  Spriggs  (15). 

These  discrepancies  in  diffusioii  coefficients  are  only  one  of  the 
factors  which  have  led  investigators  to  attempt  a  more  comprehensive  expla¬ 
nation  of  densification  phenoiriena.  A  serious  fault  of  the  Nabarro* Herring 
creep  theory  is  that  it  implies  that  grain  elongation  will  occur  during 
deformation  when  in  fact  this  is  almost  never  observed  to  occur.  The  model 
also  assumes  the  relaxation  of  shear  stresses  at  the  grain  boundary  by  a 
viscous  relaxation  phenomenon  but  the  contribution  of  this  relaxation 
mechanism  is  assumed  to  be  unimportant  in  the  deformation  process.  Creep 
studies  at  lo*.'  stress  levels  have  also  revealed  that  scratch  off¬ 
setting,  grain  rotation,  and  granular  surface  roughening  while  an  equiaxed 
grain  structure  is  maintained  can  occur  (ib).  These  observations  are  consis¬ 
tent  with  the  description  of  creep  as  a  grain  boundary  sliding  process. 

The  fact  that  a  grain  boundary  sliding  process  can  contribute  to  the 
creep  of  metals  was  reported  as  early  as  1913  (17).  However  it  was  not  until 
the  early  sixties  that  investigators  realized  that  when  creep  occurs  by  a 
diffusional  mechanism  grain  boundary  sliding  must  occur  to  maintain  specimen 
coherency  (18).  Creep  deformation  may  also  be  regarded  as  due  to  grain 
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boundary  sliding  with  diffusion  of  ions  to  accomodate  incompatabilities 
which  appear  at  the  grain  surfaces  (19).  Whichever  approach  to  grain  boun¬ 
dary  sliding  is  taken  the  two  concepts  are  physically  identical  and  when 
correctly  modeled  they  yield  precisely  the  same  result  for  the  creep  rate  (20). 

The  grain  boundary  sliding  processes  responsible  for  creep  of  poly¬ 
crystalline  solids  can  be  thought  of  in  the  following  way.  Due  to  the 
resolved  shear  stresses  acting  on  the  grain  boundaries,  adjacent  grains  may 
undergo  a  translation  with  respect  to  one  another.  However  this  translation 
(or  grain  boundary  sliding)  generates  incompatabilities  where  the  boundary 
deviates  from  a  perfect  plane.  In  order  for  sliding  to  continue  it  is 
necessary  that  an  accomodation  process  occur.  Generally  it  is  the  accomoda¬ 
tion  of  these  incompatabilities  which  controls  the  extent  and  rate  of  sliding 
(19).  The  accomodation  may  be  purely  elastic.  Elastic  stresses  build  up 
at  points  of  boundary  curvature  until  they  balance  the  applied  ones.  At 
high  temperatures  (T  ^  0.5  Tj^j)  other  accomodation  processes  are  possible. 

High  temperatures  combined  with  low  applied  stresses  are  conditions  under 
which  the  stresses  builtif>at  a  non-planar  boundary  may  be  relaxed  by  a  diffu¬ 
sive  flux  of  vacancies  from  areas  of  the  boundary  in  tension  to  those  areas 
in  compression.  And  for  conditions  of  high  temperatures  and  high  applied 
stress  levels  the  accomodation  process  may  occur  by  plastic  flow  due  to 
dislocation  motion. 

Gifkins  and  Snowden  (21)  have  examined  the  problem  of  including 
grain  boundary  sliding  contributions  to  creep  of  metals  and  metal  alloys. 

Based  on  the  idea  of  the  movement  of  double  ledges  or  protrusions  along 


the  grain  boundary  these  investigators  have  derived  an  equation  to  predict 
the  rate  of  sliding  of  a  boundary  oriented  at  45°  to  the  applied  stress. 
However  this  equation,  given  in  Table  III,  gave  sliding  rates  approximately 
30  times  larger  than  observed  experimentally  when  applied  to  the  sliding 
rates  of  lead  bicrystals. 

The  most  comprehensive  analysis  of  creep  due  to  grain  boundary 
sliding  accommodated  at  triple  points  and  grain  boundary  ledges  has  been 
presented  by  Raj  and  Ashby  (19).  These  investigators  discuss  the  process 
of  sliding  with  diffusional  accommodation  and  derive  a  general  equation  for 
the  sliding  rate  of  nonplanar  boundaries.  This  equation,  given  in  Table 
Ill,  is  based  on  the  case  where  sliding  occurs  at  a  boundary  containing 
rectangular  steps  and  where  these  steps  can  be  described  in  terms  of  their 
height,  width,  and  spacing  by  a  Fourier  series  (19). 

As  an  approximation  of  the  microstructure  of  a  polycrystalline 
solid,  a  two-dimensional,  hexagonal  array  of  grains  can  be  used.  The 
grain  boundary  shape  of  this  hexagonal  array  can  then  be  approximated  by 
a  sine  wave  cf  amplitude,  d/2  and  wavelength,  2d  where  d  is  the  grain  size 
as  shown  in  Figure  3.  By  use  of  these  simplifying  ai  ‘umptions  the  equation 
developed  by  Raj  and  Ashby  (19)  (given  in  Table  III)  can  be  used  to  obtain 
an  expression  for  the  engineering-shear  strain  rate  as  shown  in  equation  1. 
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TABLE  III 

Model  developed  by  Gifkins  and  Snowden  (21): 
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Note  that  If  the  second  term  within  the  brackets  is  <<  1  this 
equation  reduces  to  an  expression  similar  to  one  obtained  from  the  Nabarro- 
Herring  stress -directed  diffusional  model  (1,2).  Alternatively,  if  the 
second  term  within  the  brackets  is  >>1  this  equation  yields  results  similar 
to  those  obtained  by  Coble  (9)  for  grain  boundary  diffusion  controlled  creep. 

Up  until  the  present  time  it  has  not  been  possible  to  obtain  conclu¬ 
sive  evidence  as  to  the  role  which  diffusion  plays  in  the  creep  deformation 
of  ceramic  bodies.  One  of  the  major  reasons  for  this  is  the  lack  of  good 
diffusion  data  taken  on  polycrystalline  specimens.  The  use  of  an  intrinsic 
diffusion  coefficient  obtained  on  a  single  crystal  of  high  purity  by  tracer 
diffusion  techniques  is  not  an  appropriate  value  for  a  test  of  a  densifica- 
tion  model.  The  factors  known  to  influence  the  diffusion  coefficient  in 
polycrystalline  bodies  such  as  the  presence  of  grain  boundaries,  impurity 
content  and  distribution,  and  the  dislocation  density  are  not  taken  into 
account  in  making  such  a  comparison. 

In  order  to  further  our  understanding  of  densification  phenomena, 
there  is  a  need  for  a  set  of  unequivocal  diffusion  and  densification  experi¬ 
ments  designed  to  yield  diffusion  coefficients  under  experimental  conditions 
as  nearly  Identical  as  possible.  With  such  data  at  hand  it  should  then  be 
possible  to  explain  the  observed  descrepancies  between  predicted  behavior 
and  that  actually  observed  during  the  densification  process.  To  obtain  such 
information  a  program  involving  a  multicomponent  oxide,  magnesium  aluroinate 
spinel,  is  outlined  in  the  following  section. 
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Purpose  of  Investigation 

The  purpose  of  the  proposed  investigation  is  four-fold:  (1)  to 
study  the  densiflcation  kinetics  of  a  multicomponent  oxide;  (2)  to  obtain 
an  appropriate  model  for  the  densiflcation  process  and  from  the  model  to 
obtain  diffusion  coefficients  and  an  activation  energy  for  diffusion;  (3) 
to  perform  diffusion  experiments  on  the  polycrystalline  specimens  formed 
during  the  pressure  sintering  experiments  so  that  diffusion  coefficients 
can  be  obtained  on  specimens  whose  structure  and  chemical  makeup  are  as 
similar  as  possible  to  that  of  the  pressure  sintered  specimens;  and  (4) 
to  explain  in  terms  of  basic  mechanisms  the  observed  densiflcation  behavior 
of  the  multicomponent  compound. 

Method  of  Investigation 

The  compositions  to  be  used  in  this  Investigation  will  be  selected 
m  the  spinel  solid  solution  region  within  the  MgO-Al^O^  phase  diagram, 
several  features  of  this  system  make  it  attractive  for  study.  The  compound 
MgAl20^  is  regarded  as  a  model  material  and  is  representative  of  other 
spinel-type  compounds  such  as  MgCr20^,  NiAl20^,  NiCr20^,  etc.  Detailed 
studies  of  the  crystal  chemistry  of  MgA.l20^  have  been  made  and  the  locations 
of  various  cations  in  the  spinel  lattice  have  been  determined  (22,23). 

Spinel  is  of  importance  in  itself  due  to  its  applications  as  a  transparent 
armor  candidate,  as  a  high  strength-low  porosity  substrate  for  epitaxially 
deposited  silicon  integrated  circuits,  and  as  a  bonding  phase  in  high  alumina 
refractories. 


mi 
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A.  Pressure  Sintering  Experiments 

Spinel  powder  has  been  obtained  in  two  molar  ratios  of  1^0  : 

>  50:50  and  45:55.^  These  powders  will  be  used  in  the  pressure  sintering 
studies  and  some  of  the  specimens  formed  from  them  will  later  be  utilized 
in  the  diffusion  studies. 

All  pressure  sintering  runs  will  be  made  in  the  Vacuum  Industries 
hot  press  located  in  Coxe  Laboratory.  The  die  assembly  will  consist  of 
TZM  punches  and  a  die  body  having  a  graphite  liner.  Graphite  spacers  will 
be  placed  between  the  plungers  and  the  powder  so  that  the  powder  is  sur¬ 
rounded  by  a  graphite  environment. 

The  ram  displacement  during  the  run  will  be  monitored  with  an 
LVDT  driving  an  X-Y  plotter.  By  measuring  the  final  height  of  the 
specimen  the  instantaneous  strain  rate  at  any  point  in  the  run  can  be 
calculated  from  the  ram  displacement  versus  time  data  after  suitable 
corrections  are  made  for  thermal  expansion  of  the  pressure  train. 

Pressure  sintering  runs  will  be  made  at  temperatures  between  1200° 
and  1500°C  at  pressures  between  200  and  20,000  psi.  The  denslfication 

rate  will  be  determined  for  each  of  the  runs  and  the  pressure-temperature 

.  1 

regime  where  the  relationship  c  «  ka  is  obeyed  will  be  outlined.** 

The  strain  rate  data  obtained  from  this  region  will  be  analyzed  by  use  of 

*  The  powders  were  purchased  from  Dr.  Morgan  of  the  Franklin  Institute 
Research  Laboratory. 

**  Ruromler  and  Palmour  (24)  have  shown  that  the  denslfication  kinetics  of 
magnesia-rich  spinel  obeyed  this  relationship  in  the  range  1260°  to 
1450°C  and  500  to  3000  psi  applied  pressure. 
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an  appropriate  denaification  model.  A  value  of  the  apparent  diffusion 
coefficient  as  a  function  of  temperature  will  be  determined  and  the  activa** 
tion  energy  for  diffusion  will  be  calculated  from  an  Arrhenius- type 
relation. 

The  final  grain  size  of  each  specimen  will  be  measured  and  the 
strain  rate  dependence  on  grain  size  will  be  determined.  Several  pressure 
sintering  runs  will  also  be  devoted  to  a  determinaticn  of  the  grain  growth 
kinetics  of  the  specimens  so  that  if  necessary  it  would  be  possible  to 
correct  the  model  for  changes  in  grain  size  with  time. 

Some  preliminary  data  has  been  obtained  and  is  given  in  the  Appendix. 

B.  Diffusion  Experiments 

The  diffusion  experiments  are  designed  so  that  the  data  obtained  will 
compliment  that  obtained  from  the  densification  studies.  These  experiments 
will  yield  a  value  for  the  interdiffusion  coefficient  as  a  function  of  tem¬ 
perature  and  will  also  provide  information  of  the  effects  of  grain  size  on 
the  diffusion  rate, 

A  molydenum  wire  wound  tube  furnace  obtained  from  Lamont  Scientific 
Company  will  be  used  to  perform  the  diffusion  experiments.  The  oxygen  par¬ 
tial  pressure  within  the  furnace  will  be  controlled  by  passing  a  metered 
flow  of  a  CO-CO^  mixture  over  the  diffusion  couples. 

Preliminary  diffusion  anneals  are  to  be  made  with  a  single  crystal 
of  50:50*  spinel  and  a  single  crystal  of  25:75*  spinel  making  up  each 

*  These  ratios  denote  the  molar  ratio  of  MgO  :  AI2O2  in  the  costpound. 
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aandwlch'type  diffusion  couple.  These  anneals  will  be  carried  out  at 
temperatures  between  1200°  and  1700°C  for  tiires  varying  between  10  hours 
and  150  hours.  The  resultant  composition  gradient  across  the  diffusion 
couple  interface  will  be  determined  by  an  ARL  electron  beam  microprobe. 

The  data  will  be  analyzed  by  both  the  Matano  (25)  and  the  Wagner  (26) 
methods  and  the  cationic  interdiffusion  coefficient  as  a  function  of  tem¬ 
perature  will  be  calculated.  An  activation  energy  for  diffusion  will  be 
determined  by  use  of  an  Arrhenius -type  relationship. 

This  preliminary  sf'.udy  will  enable  the  investigator  to  obtain 
familiarity  with  specimen  preparation  procedures  and  techniques  of  data 
analysis.  The  best  method  by  which  to  analyze  the  data  (either  the  Matano 
or  the  Wagner  method)  will  be  determined  during  this  phase  of  the  investi¬ 
gation  and  the  results  obtained  will  be  compared  to  those  of  Whitney  (27) 
who  did  work  in  the  same  system. 

The  remainder  of  the  diffusion  experiments  will  be  directed  toward 
the  determination  of  cationic  interdiffusion  coefficients  and  activation 
energies  for  diffusion  in  polycrystalline  spinel  as  a  function  of  grain 
size  and  temperature.  Prior  to  preparation  of  the  diffusion  couples, 
pressure  sintered  specimens  of  both  compositions  will  be  annealed  at  tem¬ 
peratures  high  enough  to  cause  grain  growth  to  occur.  The  temperatures  will 
be  selected  so  that  grain  sizes  varying  by  approxiroately  two  orders  of  magni¬ 
tude  are  obtained  (for  example,  from  l^jim  to  100  (im). 

Specimens  with  a  molar  composition  ratio  of  50:50  will  be  used  as 
one  half  of  the  couple  while  the  45:55  spinel  will  be  used  as  the  other 
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half.  The  grain  sizes  of  these  two  bodies  will  be  matched  as  closely  as 
f  ssrble.  Diffusion  couples  will  be  placed  in  the  furnace  and  annealed 
at  a  temperature  ranging  from  1200°  to  1500°C.  From  analysis  of  the  com¬ 
position  gradient  through  the  couple  after  the  anneal  by  use  of  an  electron 
beam  microprobe  the  interd if fusion  coefficient  will  be  determined.  The 
specific  mathematical  analysis  to  be  used  will  be  determined  from  the 
studies  on  single  crystals  mentioned  earlier.  A  temperature-grain  size 
matrix  will  be  set  up  so  that  from  data  obtained  on  a  series  of  diffusion 
couples  having  different  g'.ain  sizes  and  having  been  annealed  at  different 
temperatures  the  effects  of  grain  size  on  the  interdiffusion  coefficient 
as  a  function  of  temperature  can  be  determined.  The  activation  energy  for 
diffusion  can  then  be  obtained  from  an  Arrhenius  relationship  between  the 
interdiffusion  coefficient  and  temperature. 

After  preliminary  data  from  the  densification  studies  and  from  the 
diffusion  experiments  are  obtained  and  analyzed  a  concittee  meeting  will  be 
called  to  discuss  the  results.  The  purpose  of  this  meeting  is  to  determine 
the  area  which  should  receive  the  most  emphasis  in  subsequent  experiments. 


APPENDIX 

Preliminary  results  have  been  obtained  from  pressure  sintering 
studies  on  45:55  (MgO  :  Al20^  molar  ratio)  spinel  at  1350°C  and  at  applied 
pressures  between  4000  and  9000  psi.  An  X-Y  recorder  with  the  Y-coordinate 
monitoring  an  LVDT  attached  to  t'.;e  hot  press  ram  and  the  X>coordinate  giving 
the  elapsed  time  was  used  to  record  the  data.  The  primary  shrinkage  data 
was  taken  from  a  dial  gauge  attached  to  the  ram  of  the  hot  press  and  whose 
smallest  division  was  0.0001  inch. 

The  strain  rate  data  obtained  from  the  densification  measurements 
was  plotted  as  a  function  of  applied  stress  at  several  porosity  levels  on 
a  log>log  scale.  As  seen  in  figure  A1  the  slope  of  the  plots  (equal  to 
the  stress  exponent,  n)  was  found  to  be  equal  to  6. 

The  high  value  for  the  stress  exponent  indicates  a  high  strain 
rate  sensitivity  on  stress.  Values  of  n  ranging  from  4  to  6  have  been 
theorized  to  indicate  that  a  dislocation  mechanism  is  operable  during 
creep.  However  it  has  been  shown  that  a  change  in  creep  mechanism  occurs 
at  lower  stresses  and  that  the  value  of  n  in  this  region  of  low  stress  is 
equal  to  1  (28). 

Secause  of  this  it  has  been  decided  that  it  is  necessary  to  keep 
the  applied  stress  below  5000  psi.  In  order  to  achieve  final  densities 
greater  than  95Z  of  theoretical  densification  studies  will  be  carried  out 
at  temperatures  between  1350“  and  1500“C. 
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ABSTRACT 

A  discussion  is  presented  of  the  stress  and  energy  cri¬ 
teria  for  brittle  fracture.  It  is  concluded  that  on  thermodynamic 
grounds  the  stress  criterion  for  catastrophic  failure  cannot  be 
satisfied  when  t^ie  energy  criterion  for  mechanical  instability 
is  met.  It  is  suggested  that  the  Griffith  criterion  defines  the 
miniminn  stress  level  for  crack  propagation  by  a  thermally  acti¬ 
vated  process.  Catastrophic  crack  propagation  can  be  achieved 
by  the  expenditure  of  unrecoverable  elastic  strain  energy  in 
addition  to  the  surface  free  energy  of  the  new  crack  surfaces. 

It  is  estimated  that  the  fracture  surface  energy  of  an  entirely 
brittle  material  is  at  least  three  times  the  surface  free  energy. 
It  is  also  concluded,  contrary  to  widespread  opinion,  that 
surface  free  energy  cannot  be  measured  by  means  of  a  fracture 
experiment. 


Introduction 

The  low  tensile  strength  of  brittle  solids  was  attributed  by 
Griffith  to  microcracks  (1).  In  the  formulation  of  his  theory,  Griffith 
hypothesizes  that  under  conditions  of  thermodynamic  reversibility,  fracture 
will  occur  when  on  crack  extension  the  change  in  energy  of  the  stress  field 
of  the  crack  equals  or  exceeds  the  surface  free  energy  of  the  resulting  new 
fracture  surfaces.  This  condition  for  crack  instability  can  be  expressed  by: 

dW/dl  >  2ys  (1) 

where  W  is  the  energy  of  the  stress  field  of  the  crack,  1  is  the  crack  length 
and  Yg  is  the  surface  free  energy  of  the  brittle  material.  Results  of 
cleavage  experiments  on  single  crystals  appear  to  establish  the  validity  of 
this  hypothesis  (2-5). 

Numerous  recent  studies,  however,  have  shown  that  in  many  brittle 
solids  even  at  low  temperature  the  energy  required  to  propagate  a  craqk 
greatly  exceeds  the  surface  free  energy,  in  apparent  contradiction  with  the 
Griffith  hypothesis.  This  discrepancy  has  been  attributed  to  surface 
roughness  and  to  energy  dissipation  processes  such  as  plastic  and  viscous 
flow  at  the  crack  tip  (6-13). 

Without  at  this  time  wishing  to  present  a  complete  review  of  brittle 
fracture,  it  is  the  purpose  of  this  coniraunication  to  suggest  that  the  work 
required  to  propagate  a  crack  in  an  entirely  brittle  solid  indeed  is  expected 
to  exceed  the  surface  free  energy.  The  excess  en^»rgy  consists  of  elastic 
strain  energy  which  must  be  provided  to  bring  the  material  to  fracture  but 
which  cannot  be  recovered  once  fracture  has  taken  place. 

Discussion  and  Theory 
Criteria  of  Fracture  of  a  Brittle  Fiber 

The  proposed  concept  of  non-recoverable  elastic  energy,  is  illvis- 
trated  most  easily  by  considering  the  criteria  for  failure  of  a  simple 
mechanical  model,  consisting  of  a  uniaxially  stressed  fiber.  The  fiber  is 
assumed  to  be  entirely  brittle  and  flaw-free  and  to  exhibit  linear  elastic 
stress-strain  behavior  up  to  the  interatomic  cohesive  fracture  stress  (o^.). 

The  load  is  applied  at  the  fiber  ends  such  that  the  stress  is  uniformly 
distributed  along  the  total  length  of  the  fiber. 
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Mechanical  stability  of  the  fiber,  under  conditions  of  thermodynamic 
reversibility,  which  forms  the  basis  for  the  Griffith  theory,  will  be  con¬ 
sidered  first.  Under  these  conditions,  at  the  critical  stress  the  total 
potential  energy  in  the  fiber  equals  the  total  surface  energy  of  the  new 
fracture  surfaces  created  on  failure.  For  failure  at  a  single  planar  cross 
section  oriented  perpendicularly  to  the  fiber  length,  this  coiidition  occurs 
when: 

L/2E  =  2Y3  (2a) 

or 

‘^cr  =  ^^‘'s 

where  is  the  critical  fracture  stress  under  conditions  of  reversibility, 

L  is  the  fiber  length  and  E  is  Young's  modulus  of  elasticity. 

Failure  of  the  fiber  can  also  be  considered  from  the  point  of  view 
of  the  energy  required  to  completely  separate  two  adjacent  atom  planes  at 
the  interatomic  cohesive  stress  (o^).  Theories  (14)  for  have  shown  that 
this  requires  a  potential  energy  per  interatomic  distance  of  the  order  of 
~Ys  P®*'  unit  area  of  fracture  plane.  At  the  critical  fracture  stress,  o^r 
of  equation  2,  it  Is  easily  calculated  that  the  potential  energy  between 
tx/o  adjacent  atom  planes  equals  2Y3  a/L,  where  a  is  the  interatomic  distance. 
Since  in  any  fracture  experiment  L>  >a,  the  energy  density  in  the  fiber  at 
a  =  is  insufficient  for  fracture  to  occur.  As  a  consequence  over  the 

stress  range: 

CT  <  cr  <  a  (3) 

cr  c 

the  fiber  is  mechanically  unstable  in  a  thermodynamic  sense,  but  fracture 
will  not  take  place  since  the  stress  condition  for  failure  is  not  satisfied. 
It  is  still  feasible,  however,  for  fracture  to  occur  over  this  stress  range 
(equation  3)  if  the  potential  eno-rgy  is  redistributed  to  the  plane  or  frac¬ 
ture  by  a  thermally  activated  process.  At  low  temperatures  and  long  fiber 
length  this  appears  unlikely. 

By  stressing  any  segment  of  the  fiber  to  the  cohesive  fracture  stress 
(c^.),  catastrophic  failure  can  be  achieved.  However,  in  view  of  the  loading 
conditions  imposed,  this  requires  stressing  the  fiber  as  a  whole  to  this 
value  of  stress  with  an  energy  expenditure,  i.e.  the  work  of  fracture  (G) 
of  the  order: 


(4) 


G  «  2y  (L/a) 
s 

On  failure,  the  fraction  of  this  energy  in  excess  of  the  surface  free 
energy  of  the  new  fracture  surfaces  is  not  recoverable  as  useful  work,  but  is 
transformed  into  vibrational  energy  (i.e.  heat,  phonons,  etc.)  in  the  now 
stress^free  fiber  segments.  It  is  absolutely  essential,  however,  that  this 
non-recoverable  elastic  strain  energy  be  provided  if  fracture  of  the  fiber 
of  the  fiber  is  to  occur. 

Criteria  of  Stability  of  a  Crack  in  a  Brittle  Solid 

Failure  by  the  propagation  of  cracks  in  a  brittle  solid  can  be  con¬ 
sidered  as  the  successive  fracture  of  atomic  bonds  in  the  flaw- free  material 
ahead  of  the  crack  tip.  Energy  and  stress  criteria  such  as  those  for  the 
previously  discussed  fiber  also  apply  to  the  stability  and  energy  expended 
in  the  propagation  of  a  crack.  Plastic  or  viscous  flov?  is  assumed  absent 
throughout  the  discussion. 

In  the  Griffith  theory  the  energy  released  on  crack  extension,  dW/dl 
of  equation  1,  is  obtained  from  the  whole  stress  field  surrounding  the  crack 
and  not  just  from  the  crack  tip  alone.  As  a  consequence,  when  the  applied 
stress  is  such  that  equation  I  is  satisfied,  under  conditions  of  thermodyna¬ 
mic  reversibility,  there  is  insufficient  energy  concentrated  at  the  crack 
tip  to  provide  the  surface  energy  required  for  crack  extension;  catastrophic 
crack  propagation  simply  cannot  occur.  On  the  basis,  then,  of  energy  consi¬ 
deration,  at  the  condition  of  equation  1  the  stress  level  at  the  crack  tip 
cannot  possibly  be  equal  to  the  interatomic  cohesive  stress.  This  conclusion 
disagrees  with  Orowan's  results  (15)  which  suggest  that  the  Griffith  energy 
criterion  also  satisfies  the  critical  stress  condition  for  brittle  fracture. 

Crack  propagation,  however,  can  occur  at  the  conditions  of  equation 
1,  if  the  energy  of  the  stress  field  of  the  crack  is  transported  to  the 
crack  tip  by  thermal  means.  This  implies  that  the  Griffith  criterion,  in 
fact  defines  the  minimum  stress  condition  required  to  cause  crack  growth  by 
a  thermally  activated  process  rather  than  the  stress  condition  for  catas¬ 
trophic  fracture. 

In  analogy  to  the  fiber,  discussed  above,  the  crack  can  be  made  to 
propagate  in  a  catastrophic  manner  even  at  a  very  low  temperature  by  raising 
the  applied  stress  to  a  value  such  that  the  stress  at  the  crack  tip  equals 
the  interatomic  cohesive  stress.  However,  as  for  the  previously  discussed 
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fiber,  the  stress  is  applied  at  a  site  far  removed  from  the  tip  of  the  crack. 
As  a  consequence.  If  the  material  at  the  cracic  tip  is  stressed  to  the  cohe¬ 
sive  stress,  material  adjacent  to  the  plane  of  crack  propagation  must  be 
stressed  as  well.  This  requires  an  expenditure  of  elastic  strain  energy  in 
addition  to  the  energy  (2Yg)  required  for  creation  of  the  new  crack  surfaces. 
It  is  the  sum  of  these  two  energies  that  constitutes  the  work  of  fracture. 

As  a  result,  the  work  of  fracture  measured  in  a  fracture  experiment  consti¬ 
tutes  the  energy  required  to  propagate  a  crack  at  a  level  of  applied  stress 
and  a  corresponding  value  of  dW/dl  (eqtiation  1)  such  that  the  stress  at  the 
crack  tip  equals  the  interatomic  cohesive  stress.  However,  under  these 
conditions  fracture  is  not  reversible  and  once  crack  propagation  occurs,  no 
mechanism  exists  by  which  the  excess  potential  energy  required  for  the 
initiation  of  crack  propagation  can  be  converted  into  useful  work.  It  is 
merely  dissipated  into  heat  or  vibration. 

It  is  of  interest  to  obtain  an  estimate  for  the  work  of  fracture  in 
terms  of  the  energy  and  stress  criteria  as  discussed  above.  For  simplicity, 
due  to  the  complex  nature  of  the  stress  distribution  near  cracks,  an  esti¬ 
mate  will  be  made  of  the  total  potential  energy  required  to  separate  two 
neighboring  atom  planes  at  the  interatanic  cohesive  stress  with  a  stress 
field  adjacent  to  the  plane  of  fracture  similar  to  the  distribution  of  stress 
near  the  tip  of  a  crack.  As  suggested  by  the  Airy  stress  function  (16), 
within  a  distance  approximately  one  tenth  of  the  crack  size,  the  stress  field 
near  the  crack  tip  can  be  expressed  by: 

o  =  ky"^  (5) 

where  k  is  a  constant  and  y  is  the  distance  from,  a  plane  midway  between  the 

atom  planes  being  separated.  The  constant  k  can  be  evaluated  from  the  boun- 

+k 

dary  condition:  when  y  =  a/2,  o  =  a^..  This  gives  k  =  (a/2)  ^  and: 

CT  =  (a/2y)^  (6) 

Taking  the  stress  distribution  symmetric  about  y  =  0  and  noting  that  for 

-a/2<y<a/2  the  work  done  equals  2ys>  the  total  work  in  separating  unit 
area  of  atom  plane  (i.e.  the  work  of  fracture)  becomes: 

R 

G  =  2Y3+  5  (Oj,'  a/2yE)  dy  (7) 

a/2 
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where  R  represents  the  "range  of  influence",  of  the  stress  field  at  the  tip 
of  the  crack.  Setting  =  (Vg  E/a)^  after  (14),  followed  by  integration 
yields: 

G  =  2Yg+%Y3'f'n  (2R/a)  (8) 

To  obtain  a  value  for  G  that  is  more  representative  of  the  stress  field  of 
a  crack,  the  multiaxlal  nature  of  the  stress  field  can  also  be  considered. 
Cook  and  Gordon  (17)  showed  that,  but  for  the  immediate  vicinity  of  the  tip 
of  a  crack  in  a  flat  plate  the  stress  is  of  a  uniform  biaxial  nature.  This 
introduces  an  additional  factor  2(l-v)  in  the  second  terms  of  equation  8, 


where  v  is  Poisson's  ratio,  such  that: 

G  =  2y„  +  (1-v)  Y,  't.n(2R/a)  (9) 

s  s 

An  estimate  of  G  may  be  made  on  the  basis  of  a  numerical  example.  Conserva** 
tlvely  estimating  the  range  of  influence,  R,  to  extend  over  only  one 
hundred  Interatomic  distances  (i.e.  R  =  100a),  gives  for  v  =  0.25; 

G  «  6.0  Y„  (10) 

s 

or,  if  preferred  in  terms  of  the  fracture  surface  energy,  Yj  ®  G/2; 

V,  «  3.0  (11) 


This  result  suggests  that  in  propagating  a  crack  the  additional 
energy  required  to  stress  material  located  within  only  one  hundred  inter¬ 
atomic  distances  already  has  the  effect  of  nearly  tripling  the  fracture 
surface  energy.  Since  only  one  hundred  interatomic  distances  were  consi¬ 
dered,  equation  11  clearly  represents  an  underestimate  for  Yf*  addition, 
in  polycrystalline  materials  surface  roughness  can  easily  lead  to  a  further 
doubling  or  tripling  of  Yf  to  a  value  perhaps  as  high  as  10  Regardless 

of  the  actual  value,  it  is  clear  that  the  non-recoverable  elastic  strain 
energy  can  cause  the  fracture  surface  energy  in  an  entirely  brittle  material 
to  exceed  the  surface  free  energy  by  an  appreciable  amount.  As  a  conse¬ 
quence,  no  need  exists  to  resort  to  energy  dissipation  by  plastic  or  viscous 
flow  to  completely  explain  the  observed  discrepancies  between  surface  free 
energy  and  fracture  surface  energy.  Clearly,  energy  dissipative  processes 
such  as  plastic  or  viscous  flow  can  contribute  further  to  the  work  of  frac¬ 
ture  but  they  are  not  the  sole  mechanisms. 

The  present  discussion  leads  to  the  Inevitable  conclusion  that  sur¬ 
face  free  energy  of  a  brittle  solid  cannot  be  measured  by  means  of  a 
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fracture  experiment  using  specimens  of  macroscopic  size.  Similar  conclusions 
were  also  expressed  by  Bikerman  (18)  and  Thompson  et  al.  (19).  Those  studies 
in  which  fracture  experiments  result  in  apparent  values  of  surface  free 
energy  need  to  be  re-examined.  Surface  free  energy  in  principle,  can  be 
measured  by  mechanical  means  by  determining  the  minimum  stress  required  for 
crack  growth  by  a  thermally  activated  process.  At  low  temperature  such 
experiments  are  expected  to  be  time  consuming. 

The  surface  free  energy  (Vg)  is  directly  proportional  to  Young's 
modulus  of  elasticity  as  shown  by  Gilman  (20).  As  a  consequence  from  equa¬ 
tion  11,  which  relates  Vf  ^nd  Ys»  fracture  surface  energy  is  also 
expected  to  be  proportional  to  Young's  modulus.  Observations  by  Wlederhom 
(10)  on  the  fracture  surface  energy  of  various  glasses  with  different  values 
of  Young's  moduli  are  in  support  of  this  hjrpothesis. 

An  additional  observation  must  also  be  made  in  relation  to  the 
interpretation  of  brittle  fracture  data.  Since  equation  1  represents  the 
minimum  energy  condition  required  to  render  a  crack  thermodynamically 
unstable  and  equation  8  or  11  gives  the  energy  required  for  catastrophic 
crack  propagation,  over  the  range: 

2y  <  dW/dl  <  G  (12) 

s 

the  crack  is  unstable  in  a  thermod3mamic  sense;  catastrophic  fracture  will 
not  occur  but  crack  growth  can  take  place  by  means  of  a  thermally  activated 
process.  As  a  result,  in  any  fracture  experiment  carried  out  at  a  tempera¬ 
ture  in  excess  of  0®K,  some  subcritical  crack  growth  prior  to  catastrophic 
failure  is  to  be  expected. 

The  thoughts  expressed  in  the  present  note  can  be  used  to  qualita¬ 
tively  interpret  the  observations  of  Congleton,  Fetch  and  Shiels  (21)  for 
the  temperature  dependence  of  the  fracture  surface  energy  of  polycrystalline 
aluminum  oxide.  At  -196®C,  these  authors  observed  Yf  ^  3  x  10^  ergs. 
cm"2,  decreasing  to  approximately  1.5  x  10^  ergs.  cra“2  at  250‘’C  followed 
by  a  rapid  increase  to  about  6  x  10^  ergs.  cm"2  at  500®C.  Taking  Ys  ^ 

2,000  ergs.  cm”2  and  assuming  that  surface  roughness  is  such  that  the  actual 
surface  area  is  doubled,  equation  11  suggests  a  Yf  ^  lx  10^  ergs.  cm"2 
or  even  higher.  As  a  result,  at  -196'’C  perhaps  as  much  as  one  half  of  the 
value  of  Yf  be  attributed  to  non-recoverable  elastic  energy.  The 
observed  decrease  in  Yf  from  -lOb^C  to  250’’C  can  be  attributed  to  the 


redistribution  of  elastic  energy  from  the  stress  field  of  the  crack  to  the 
crack  tip.  This  reduces  the  applied  stress  required  to  raise  the  crack  tip 
stress  to  the  cohesive  strength  which  manifests  itself  in  an  apparent 
decrease  in  the  fracture  surface  energy,  (kily  the  increase  in  Vf  at  temper* 
atures  above  250‘’C,  in  the  opinion  of  the  present  writers,  is  due  to  plastic 
work  as  the  result  of  increased  dislocation  mobility.  Similar  observations 
for  the  temperature  dependence  of  the  fracture  surface  energy  have  been 
reported  by  Coppola  and  Bradt  for  silicon  carbide  (22)  and  silicon  nitride 
(23)  and  can  possibly  be  explained  by  the  above  argument. 

Summary 

The  present  discussion  points  out  that  at  the  minimum  energy  require¬ 
ment  for  mechanical  instability  under  conditions  of  thermodynamic  reversibi¬ 
lity,  stress  requirements  for  catastrophic  failure  cannot  be  satisfied. 
Catastrophic  fracture  of  an  entirely  brittle  solid  requires  an  expenditure 
of  energy  well  in  excess  of  the  surface  free  energy,  which  is  not  recoverable 


on  fracture. 
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The  Griffith  Approach  to  Brittle 
Fracture  and  Plastic  Flow  in  Glass 


Abstract 


The  Griffith  approach  to  brittle  fracture  is  considered 
with  the  emphasis  on  fractuie  surface  energy.  The  reversible 
concept  of  brittle  fracture  is  modified  to  include  the  ir re vers 
ible  processes  accompanying  the  creation  of  new  fracture 
surfaces.  These  processes  are  reviewed  with  regard  to  their 
effect  on  the  true  surface  fracture  energy.  A  model  for 
"Dislocation-like”  behavior  is  proposed  and  independent  ex¬ 
perimental  support  of  this  proposal  is  discussed.  \ 
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Introduction 


Glass  is  an  isotropic  material  with  an  amorphous  structure. 
It  breaks  in  a  brittle  way  with  little  or  no  apparent  plastic 
deformation  taking  place.  Since  glass  is  observed  to  fracture 
under  a  critical  uniaxial  tensile  stress,  the  maximum  tensile 
stress  is  the  criterion  for  failure  The  fracture  is  always 

observed  to  be  initiated  from  the  surface  for  glass  free  of  in¬ 
clusions.  The  surface  control  of  the  fracture  of  glass  is  not 
surprising  since  sources  of  internal  weakness  associated  with 
grain  boundaries  are  not  present. 

Theoretically,  the  fracture  strength  of  glass  should  equal 
the  cohesive  strength  for  breaking  bonds (i.e.,  E/10  where 
E  is  Young’s  modulus).  This  means  that  for  most  glasses,  the 

theoretical  fracture  strength  should  be  on  the  order  of  750,000- 

f3  4) 

1,500,000  psi.  Several  investigators^  *  '  have  observed 
strengths  of  this  magnitude.  However,  the  observed  strengths 
of  most  glasses  are  reported  to  be  between  10,000  and  50,000  psi. 
Griffith^^’^^  proposes  that  the  difference  between  theoretical 
strengths  and  those  experimentally  observed  is  due  to  the 
presence  of  small  cracks  which  act  as  stress  concentrators. 

When  the  material  is  stressed,  the  theoretical  strengfth  is 
reached  in  small  voltones  of  glass  even  though  the  average  stress 
is  low. 
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Griff i til  Approach 

Griffith^^^  proposes  that  the  unstable  propagation  of  a 
crack  must  result  in  a  decrease  in  free  energy  of  the  system. 
The  source  of  energy  for  crack  propagation  is  the  release  of 
strain  energy  which  goes  into  creating  the  new  surface.  The 
total  free  energy  of  the  system  containing  a  crack  is  of 
the  form: 


=  M  -  Pg  + 


P  = 


2„2 

p^  =  II£C 


Pg  =  4Cys  (4) 

where;  p  =  elastic  strain  energy  when  no  cracks  are  present  (p5^f(c)  ) 

p  =  stored  elastic  strain  energy  which  is  released  as 
the  crack  grows 

Pg  =  total  surface  energy  resulting  from  the  crack 
a  =  applied  stress 
C  =  crack  length 

Y  =  true  surface  free  energy  (work  done  in  forming  unit  area 

^  of  nev  surface.^s  a  result  c-f  breaking  atomic bdhdS) 

V  =  volume 


The  condition  for  crack  growth  is; 


-  Ug)  so  (5) 

This  equation  indicates  that  a  crack  will  grow  when  the  incremental 
release  of  elastic  strain  energy  is  greater  than  the  incremental  in¬ 
crease  of  surface  energy  produced  as  the  crack  surfaces  are  gene¬ 
rated.  From  equation  4  the  Griffith  equation  is  obtained 

o  fracture  =  o  r 


Ill  d  moiilfiod  form,  y^  replaced  by  y,  the  fracture  surface 

o 

energy,  (energy  of  formation  of  a  unit  area  created  by  the 

fracture  process^^^).  The  reason  for  this  substitution  will 

become  evident  in  the  following  discussion. 

(3  5  ) 

Griffith  ’  considers  that  fracture  is  a  reversible 
process.  He  indicates  that  the  strains  must  be  elastic  and 
that  cracks  can  be  propagated  or  healed  by  increasing  or  de¬ 
creasing  the  load.  The  energy  absorbed  during  the  surface 
formation  is  lost  upon  surface  annihilation  (healing).  However, 
fracture  is  usually  accompanied  by  irreversible  processes  such 
as  plastic  flow,  viscous  flow  and  irreversible  chemical  reac¬ 
tions.  Wiederhorn^^^  indicates  that  these  irreversible  process¬ 
es  cause  the  fracture  surface  energy  to  be  more  or  less  than 
the  true  surf ice  free  energy.  Plastic  and  viscous  flow  cause  a 
shear  displacement  of  material  which  prevents  rehealing.  Shear 
displacements  increase  the  fracture  surface  energy  because  energy 
is  expended  during  the  deformation  process.  Stress  corrosion 
prohibits  rehealing  by  removing  material  and  chemically  altering 
the  now  surf acc-^-  -Stress  corrosion*  decreases  the  fracture  surface 
energy  because  the  chemical  reaction  supplies  energy  for  surface 
formation . 

1  H. 

From  the  relation  a  (C)^  j  *  Griffith^^^  using  e>:- 

u 

perimental  values  determined  that  the  average  value  of  a(C) 
h  /2Eys\^  % 

was  239  psi  (in)  andf — was  134  psi  (in)  .  The  discrepancy 
can  be  accounted  for  by  considering  that  the  true  surface  free 
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energy  does  not  take  irreversible  processes  into  consideration. 

The  value  of  Vg  used  by  Griffith  was  0.003  Ib/in.  Wiederhorn^^^ 

measured  the  fracture  surface  energy  of  soda-lime-silica  glass 

(similar  to  the  glass  that  Griffith  used)  to  be  about  0.021  Ib/in 
2 

(3.82  joules/m  )  which  reflects  the  irreversible  processes. 
Griffith’s  value  of  surface  energy  (extrapolated  to  room  tempera¬ 
ture  from  high  temperature  surface  tension  data)  is  only  about 
15%  of  the  measured  value.  In  addition  to  not  including  irrevers¬ 
ible  processes,  another  possible  large  error  is  suggested  by 
Morey^^*^\  At  higher  temperatures  the  structural  elements  of 
the  glass  are  mobile  and  an  equilibrium  structure  is  present 
which  has  a  minimum  surface  free  energy.  At  lower  temperatures, 
the  glass  components  are  not  mobile  and  the  nonequilibrium  struc¬ 
ture  is  characterized  by  a  higher  surface  free  energy  than  ex¬ 
trapolation  would  yield. 
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Plastic  Deformation 

Recently  reported^^^  values  of  glass  strength  are  highly 
reproducible  and  probably  indicate  a  pristine  surface  condition. 
However,  these  values  are  still  well  below  the  theoretical  co¬ 
hesive  strength.  Also  it  has  been  observed  that  calculated 
fracture  surface  energies  are  as  much  as  50  times  greater  than 
the  surface  free  energy^®  \  Marsh^^^  indicates  that  an  elastic- 
plastic  theory  could  account  for  the  discrepancies.  This  theory 
can  be  illustrated  by  considering  that  the  random  structure  of 
glass  causes  it  to  be  a  non-work  hardening  solid.  Glass  has  no 
stabilizing  mechanism  at  the  yield  point  and  so  must  fail  cata¬ 
strophically  when  the  applied  stress  equals  the  flow  stress.  The 
apparent  complete  brittleness  of  glass  is  a  result  of  the  fact 
that  fracture  occurs  at  the  flow  stress  with  negligible  elongation. 
Marsh  demonstrates  plastic  flow  around  indentations  and  scratch 
marks  on  glass  specimens.  He  indicates  that  the  plastic  flow 
around  indentations  results  from  very  high  local  stresses.  Similar 
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ii.igh..stnesses..ey.ist— iP  the  vicinity  of  *crdck  tip  just  priori to 


failure  and  therefore  plastic  flow  should  occur. 

Marsh^^^  proposes  that  the  fracture  criterion  is  that  of 
a  critical  plastic  zone  size.  Wiederhorn^^^  estimates  the 
plastic  zone  size  using  Dugdal*s  model  of  plastic  flow  at  the 
crack  tip  (see  Figure  1). 


The  following  assumptions  are  made  in  the  model^  first 
the  plastic  zone  extends  a  length  R  in  front  of  the  crack  tip; 
second,  the  material  within  the  plastic  zone  is  at  the  yield  stress 
of  the  solid.  The  plastic  zone  is  given  by' 

R  =  (7t/8)  (7) 

where  is  the  critical  stress  intensity  factor  and  Oyg  is 
the  yield  stress.  The  displacement  at  the  crack  tip  is^^^: 

V(c)  =  4  R/ttE  (8) 

Energy  absorption  is  equal  to  aygV(c).  Wiederhorn^^^  estimates 

the  size  of  the  plastic  zone  for  soda- lime-silica  glass  to  be 

-9  -10 

2.9  X  10  m  and  the  crack  opening  displacement  to  be  4.5  x  10  m. 

He  indicates  that  the  dimensions  are  not  much  larger  than  the 

silicon-oxygen  distance  of  1.6  x  10  ^^m  or  the  oxygen-oxygen 

distance  of  2.6  x  10  The  very  brittle  nature  of  glass  is 

a  consequence  of  the  small  plastic  zone  size.  Wiederhorn^^^ 

found  that  the  energy  absorption  for  a  moving  crack  is  about 

4.5  joules/m  for  soda-lime-silica  glass  as  compared  to  about 
4  /  2 

.ID  .  jiMiles/m.  for  meta.'ls  for  vhir.h-the--5‘l5Scic  zonp  is  A.uch -larger-  - 

than  for  glass.  As  a  result  for  glass,  as  a  crack  begins  to  move 

it  tends  to  propagate  catastrophically.  The  energy  absorption 
o 

of  4.5  joules/m*"  is  approximately  equal  to  the  measured  surface 
fracture  energy.  Although  this  might  seem  conclusive  for  the 
plastic  flow  argument,  Wiederhorn^^ ^  points  out  that  the  model 
uses  a  continuum  mechanics  approach  and  its  validity  is  in  question 
for  such  small  plastic  zone  sizes. 
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In  recent  work,  Wiederhorn^^^  shows  that  cracks  in  soda-lime- 
silica  glass  can  heal.  Cracks  can  heal  only  if  proper  alignment 
can  be  obtained  when  the  crack  faces  are  rejoined-  The  healing 
process  is  caused  by  atomic  forces  pulling  the  surfaces  together. 
Distortions  from  plastic  flow  will  keep  the  surfaces  apart.  There¬ 
fore,  for  cracks  to  heal,  the  crack  tip  displacement  must  be  within 
the  range  of  atomic  forces  in  glass,  i.e.,  less  than  10x10 
Considering  the  Dugdale  model,  the  crack  tip  displacement  is  about 
4.5x10  ^®m  which  is  within  the  range;  therefore  plastic  deformation 
can  occur  and  not  prevent  the  healing  process. 

Several  other  theories  have  been  proposed  to  explain  the  flow 
of  glass  around  indenters  and  account  for  the  high  fracture  surface 
energy.  Viscous  flcn-;  is  a  possible  mechanism^ The  strain 
field  caused  by  the  high  stresses  at  the  crack  tip  can  cause  a 
localized  heating  thus  allowing  the  glass  to  flow  viscously. 

Lange  has  recently  calciilated  the  degree  of  heating  at  the 
crack  tip  and  found  it  to  be  about  200°C.  This  is  well  below  the 

region  where  vi'scou's* f iov?  tfouicl'be  a  I'easonable “factor. 
fl2j 

Roughness'  ■'  of  the  surface  must  also  be  considered  to  explain 
the  h-igh  fracture  surface  energy.  If  the  fracture  surface  is 
rough,  shows  forking,  or  sub-surface  cracks  are  formed,  these  will 
contribute  to  the  fracture  surface  energy.  The  Griffith  approach 
has  a  thermodynamic  basis  and  irreversible  losses  to  the  surround¬ 
ings  owing  to  the  second  law  of  thermodynamics  accompany  the 
fracture  process 


suggests  a  qualitative  model  to  explain 
the  Insults  of  indentation  hardness  tests.  The  application  of 
hydrostatic  pressure  to  a  glass  results  in  an  elastic  compression 
and  the  recovery  is  complete  when  the  load  is  removed.  If  while 
the  glass  is  under  compression,  it  is  also  subjected  to  a  shear, 
interlocking  of  the  network  will  result  as  indicated  in  Figure  2 . 
This  entanglement  causes  an  apparent  densification  when  the  pres¬ 
sure  is  released. 

The  geometry  of  a  pyramidal  indentor  gives  rise  to  shear 
stresses  when  it  is  penetrating  the  sample .  Thj  s  suggests  that 
densification  by  high  pressure  may  be  involved  in  the  observed 
flow  resulting  from  the  indentation  tests.  Mackenzie indi¬ 
cates  that  partial  recovery  of  the  indentation  volume  flow  is 
observed  at  temperatures  much  lower  than  the  glass  transition 
temperature  (temperature  where  the  glass  can  viscously  relax). 

He  indicates  that  the  recovery  is  difficult  to  explain  if  it  is 
assumed  that  the  flow  is  caused  by  a  viscous  or  a  plastic  flow 
mechanism.  Mackenzie  also  observes  that  densifiealsion  due  to 
high  pressure  recovers  at  temperatures  much  below  the  glass 
transition.  However,  the  fact  that  only  partial  recovery  is  ob¬ 
served  indicates  that  perhaps  plastic  flow  is  also  occurring. 
Ernsberger^^^^  using  optical  techniques  has  confirmed  that  densi¬ 
fication  does  take  place  at  the  indentation.  The  experimental 
results  do  not  rtile  out  the  possibility  that  densification  and 
plastic  flow  can  occur  together.  It  should  be  pointed  out  that 
for  plastic  or  viscous  flow,  the  volux.ie  is  conserved  and  for 
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shows  experiment- 


densification  the  volume  disappears.  Taylor^^^^ 
ally  that  the  material  displaced  by  the  indentor  which  forms  the 
raised  lip  is  equal  to  the  volume  of  the  indentation.  This  conser¬ 
vation  of  volume  supports  a  viscous  or  plastic  flow  mechanism.  In 
view  of  the  temperature  limitations  calculated  by  Lange,  plastic 
flow  is  the  more  probable  mechanism. 

After  consideration  of  all  the  proposed  theories  of  flow, 
the  experimental  and  theoretical  evidence  indicates  that  plastic 
flow  is  a  possible  energy  absorbing  mechanism  in  glass  fracture . 

Discussion  of  the  mechanism  of  plastic  flow  is  complicated 

C19) 

by  the  absence  of  dislocations  in  glass.  However,  Levengood 
observes  flaw  patterns  in  glass  which  behave  under  the  influence 
of  stress  and  etching  solutions  in  a  manner  similar  to  disloca¬ 
tions  in  crystals.  also  observes  that  etch  pits  in  the 

glass  under  an  applied  stress  act  as  a  source  for  "dislocation¬ 
like"  loops.  From  the  dimensions  of  the  source,  he  calculates  a 

f  21) 

Burger’s  Vector.  Levengood  calculates  that  the  energy  to 
form  a  flaw  is  small  *ciiid  tH«C  they  car  easily  be*rormeB“at  stress 
levels  much  below  that  of  fracture.  By  a  plastic  flow  mechanism 
these  flaws  can  combine  to  form  a  microcrack  and  lead  to  subse- 
quent  fracture.  March'  '  and  Ainsworth'  ’  suggest  a  "disloca¬ 
tion-like"  mechanism.  Ainsworth  shows  that  the  flow  of  glass  is 
related  to  the  proportion  of  network  modifiers  in  the  structure. 

The  addition  of  network  modifiers  introduces  non-bridging  oxygens 
which  provide  "loose  ends"  in  the  structure  and  allows  flow  by  a 
mechanism  analogous  to  half  planes  which  terminate  at  dislocations 

^9- 
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in  crystals. 

An  attempt  to  expand  on  the  dislocation- like  mechanism  is 

given  in  the  following  discussion.  Figure  3  shows  the  Zachariasen 

typo  structure  for  silica  glass.  Also  shown  is  the  effect  of 

vulding  a  network  modifier,  in  this  case,  sodium.  The  modifier 

hiHMks  the  c.ontimious  linking  of  oxygen  tetrahedra  and  weakens 

Liu;  gla:;:;.  I’or  tlislocation  movement  in  a  crystal,  the  extea  half 

plane  i;f  atoms  sits  at  an  energy  maximum.  Forces  to  move  the 

dislocation  in  one  direction  are  balanced  by  forces  in  the  other. 

Therefore,  the  force  required  to  move  the  dislocation  is  small 

f  25  ^ 

and  plastic  flow  is  readily  observed^  A  similar  argument  can 
be  applied  to  the  modified  glass  structure  pictured.  The  non¬ 
bridging  oxygens  have  a  negative  charge  while  the  sodium  ions  have 
a  positive  charge.  As  a  rough  approximation,  the  attractive  forces 
exei’ted  by  Na^  labeled  1  will  be  balanced  by  forces  exerted  by  Na^ 
labeled  2.  However,  unlike  the  crystal,  in  which  the  dislocations 
move  by  glide  on  a  slip  plane,  to  move  the  oxygen  tatrahedra  re- 
quires  a  great  deal  of  force  to  distort  the  structure  and  several 
bonds  must  be  broken.  (For  crystals,  unit  movement  requires  only 
that  one  bond  be  broken) .  The  movement  of  non- bridging  oxygens  is 
one  of  tearing  which  follows  no  particular  plane.  The  dislocation 
line  is  pictured  in  Figure  3  and  will  have  a  variable  Burger  Vector 
as  discussed  by  Gilman^ Since  the  flow  stress  is  the  stress 
at  which  numerous  bonds  are  broken  by  a  tearing  mechanism,  then 
it  is  also  expected  that  this  will  be  the  fracture  stress.  This 
is  in  agreement  with  March* s^^^  work.  As  an  example,  Marsh  shows 
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che  flow  stress  and  the  fracture  stress  for  soda-lime-silica  glass 
are  510,000  psi,  and  520,000  psi  respectively.  From  the  above  argu¬ 
ments,  in  the  absence  of  modifiers  in  the  structure,  plastic  flow 

would  not  occur.  Then  fused  silica  (with  no  modifiers)  would  be  | 

! 

expected  to  fractui-e  in  a  manner  without  plastic  flow.  This  con-  ■! 

(  21 )  • 
elusion  is  consistant  with  Levengood’s  work  in  which  he  demon- 

strates  that  fused  silica  is  more  brittle  than  glasses  containing  ! 

(27') 

modifiers .  Schardin'-  ■'  has  also  demonstrated  this  point .  He 

found  that  the  ratio  of  the  velocity  of  brittle  cracks  to  the  shear  ; 

wave  velocity  decreased  as  network  modifiers  were  added  to  the  glass 

indicating  a  dissipative  energy.  For  fused  silica  the  ratio  was 

the  highest  revealing  the  extreme  brittleness. 

Summary ; 

The  Griffith  approach  to  brittle  fracture  explains  the  differ¬ 
ence  between  theoretical  cohesive  strength  and  observed  fracture 
strength.  Developments  in  the  field  of  fracture  mechanics  have  helped 

to  understand  the  discrepancies  between  fracture  surface  energy  and 

•  •  •  •  .  *  «  , 

surface  free  energy.  The  original  concept  of  fracture  being  a  re¬ 
versible  process  must  be  modified  to  include  the  irreversible  pro¬ 
cesses  which  several  investigators  show  to  exist.  For  glass 
fracture,  plastic  flow  should  be  considered  and  the  thought  that 
glass  breaks  in  a  completely  brittle  way  is  certainly  in  doubt.  It 
seems  unlikely  that  an  entire  technology  will  develop  around  dislo¬ 
cation-like  movement  in  amorphous  solids  as  has  happened  for 
crystalline  solids.  However,  it  is  interesting  to  speculate  how 
flow  can  occur  by  dislocation-like  mechanisms. 
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Figure  1.  The  Dugdale  model  of  elastic- plastic  deformation 
near  a  crack.  (After  Wiederhorn  C^)  ) 


Figure  2.  Mackenzie's  model  of  network  entanglement  and 
apparent  densification  under  the  influence  of 
an  applied  shear  stress.  (After  MacKenzie^^^'  ) 


Figure  3.  Zachariasen  model  for  the  structure  of  glass. 

A.  Fused  Silica  (After  Zachariasen  (28)  ) 

B.  Soda-Silica  Glass  (After  Warren  (24)  ) 
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